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' Abstract. We present the results of a spectroscopic analysis of the Trapezium cluster stars inside the Orion 

Nebula. The rotational velocities have been obtained using Fourier analysis method, finding agreement with 
. values derived from the usual method, based on linewidth measurements. The rotational velocity derived for 9^ 

' Ori C through this method is consistent with the variability of some of its spectral features that have a period of 

fH I 15.42 days. By means of the fit of H , He I and Hell observed profiles with Fastwind synthetic profiles, stellar 

parameters and wind characteristics have been derived. This methodology let us estimate the errors associated 

1 ' with these parameters. It is found that macroturbulence effects have to be included for a good fit to the Hel-ll 
O ^ lines in the spectrum of 6^ Ori C. 

By means of a very accurate study, oxygen abundances have been derived for the three B0.5V stars 9^ Ori A, 
^ ' D and 9^ Ori B. Final abundances are consistent with the nebular gas-phase results presented in Esteban et al. 

H2UU4II and are lower than those given by Cunha & Lambert II1994II . Our results suggest a lower dust depletion 
factor of oxygen than previous estimations for the Orion nebula. 

' Key words. Stars: abundances - Stars: early-type - Stars: fundamental parameters - Stars: atmospheres - Stars: 

individual: 9^ Ori C - Stars: rotation - ISM: abundances - ISM: H II regions - ISM: individual: Orion nebula 



1. Introduction 

New developments of massive star model atmosphere 
codes have led to interesting new possibilities for stellar 
spectroscopic studies. Improvements in computational 
methods as well as an increase of the efficiency of com- 
puters, have made possible the modeling of atmospheres 
of hot luminous stars, taking into account not only 
strong NLTE effects and hundreds thousands metallic 
lines producing the so-called line-blanketing (Hubeny & 
Lanz I1995|l . but also winds with expanding spherical 
geometries (Santolaya-Rey et al. 119971 Hillier & Miller 

Pauldrach et aL OnUTl Puis et al. 12003) • 
The new improvements that have been included in these 
latest generation models call for a review of previous 

Send offprint requests to: ssimon@iac.es 

* The INT is operated on the island of La Palma by the RGO 
in the Spanish Observatorium of El Roque de los Muchachos 
of the Instituto de AstroHsica de Canarias. 



results. For example, the papers by Herrero et al. H2002fl . 
Crowther et al. H2002|l . and Martins et al. (|^D7) showed 
that the SpT - Tes calibrations used previously (Vacca 
et al. 11996(1 needed to be revised to lower effective 
temperatures for a given spectral type. Recent analyses 
by Repolust et al. I|2004|l and Martins et al. H2005i in the 
Milky Way, and by Massey et al. ( 200^ at SMC 

and LMC metalicities, reinforce this result and together 
implies a need to revisit the ionizing flux distribution 
that is used for the study of H ii regions and starbursts. 

Recent work by Trundle et al. H2002|l and Urbaneja 
et al. H2005(l have shown that abundance gradients in 
some spiral galaxies derived from stellar and nebular 
studies tend to be coherent but very dependent on the 
calibration used in the strong line nebular methods. 
However, until now, there have been no detailed sys- 
tematic studies which compare results from nebular and 
stellar studies. This is the first of a series of papers 
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Fig. 1. Atlas of the int+ids spectra in the 4000-5000 A region used for this study, r Sco spectrum is from CASPEC. 
Several narrow nebular emission lines appear in the spectra of the Orion stars. In these stars, the cores of the hydrogen 
Balmer lines are contaminated by the nebular lines. Nebular and [Oiii] AA 4960,5007 lines have been arbitrarily 
diminished in the 9^ Ori A and D plotted spectra for clarity reasons. Many metal lines can be distinguished together 
with the broad H and He lines in the 09 - BO. 5 spectra. The high rotational velocity of 9^ Ori A makes the metal 
lines of this 09V star shallower, broader and blended. The known inverted P-Cygni He ii A4686 profile can be seen in 
the 6^ Ori C spectrum. 



aimed at this subject. We have selected some galactic Hii 
regions for a detailed study of the interaction between 
massive stars and the surrounding ISM, looking for 
consistency of derived parameters ( Tcs , luminosities and 
ionising flux distribution of the stellar population) as 
well as abundances of C, N, O, Si and Mg. For this first 
study we have selected the Orion nebula, a well studied 
and resolved Galactic Hii region powered by a cluster of 
a few massive stars, the Trapezium cluster. 

The Orion complex contains the massive on-going 
star forming region closest to Earth, at about only 450 
pc. The Orion nebula, M42, is part of this complex. It 
is a well known H ii region (e.g. O' Dell t2001 Ferland 
1200111 ionised by the Trapezium cluster stars {9^ Ori), 
a group of early type stars located in the core of the 
nebula. Together with 9^ OriC (HD 37022, 07V), the 
main ionising source, we find other B0.5V stars, perfect 



targets for a stellar abundance study. 
The most recent study of the chemical composition of 
the Orion nebula has been presented by Esteban et al. 
()20()4|l who reanalysed the emission line spectrum of 



the nebula to determine the physical conditions and 
abundances of the ionised gas-phase. Cunha & Lambert 
II19921 11994(1 included some of the Trapezium cluster stars 
in a survey of B-type stars in the Orion OBI association. 
They presented a spectroscopic analysis of these stars for 
determining C, N, O, Si and Fe stellar abundances. 

This paper is focused on the spectroscopic analysis 
of the Trapezium cluster stars for deriving their stellar 
parameters and oxygen abundances. The stellar parame- 
ters obtained for 9^ Ori C will be used in future papers 
as input for the modeling of the Orion Nebula with 
photoionization codes. The derived stellar abundances 
are compared to those obtained by Esteban et al. 
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through nebular studies. 

Our paper is structured as foUows: In Sect. |21 we 
present the observations. In Sect. 13 and ^ we obtain the 
wsini and the stellar parameters of our targets. The 
study of 0^ Ori C is presented in Sect. El and the oxygen 
abundance analyses in Sect. El A discussion of the results 
and the conclusions of this work are presented in Sect. [3 
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Fig. 2. Atlas of the int+ids spectra in the Hq region used 
for this study, r Sco spectrum is from CASPEC. Nebular 
Ha and [Nil] emission lines have been diminished in the 
6^ Ori A and D plotted spectra for clarity reasons. Two 
emission components can be distinguished in the 9^ Ori 
C Hq line: the narrow one is associated with the nebula, 
the other one is variable and it is associated with the star 
(see sectionEJ. The blue wing of the Hq line is affected by 
some cosmetic in the iNT observations. 



2. The observations 

The bulk of the observations used here were carried out 
with the Isaac Newton 2.5m. Telescope (int) at the 



Roque de los Muchachos observatory in La Palma on 
the 21th of December 2002. The Intermediate Dispersion 
Spectrograph (ids) was used with the 235 mm camera 
and two different gratings. We have observed the spectral 
region between 4000 and 5050 A using the H2400B 
grating, which resulted in an effective spectral resolution 
R~7500 (equivalent to a 0.23 A/pixel resolution and ^ 
2.6 pixel FWHM arc hues). For the Hq region the H1800V 
grating was used, resulting in a similar spectral resolution 
(0.3 A/pixel, R~8000). With the configurations used 
three exposures were needed to cover the whole range. 
A large number of flat fields and arcs for the reduction 
process were obtained. 

The reduction and normalization of the spectra was 
made following standard techniques, with IRAF and our 
own software developed in idl. The signal-to-noise ratio 
(SNR) of the reduced spectra depends on the spectral 
range, but is usually about 200-250 in the blue region 
and 250 in the region (see Table |2l. 

We found some problem when rectifying the int 
observations in the region. In the blue wing of the 
Hq line we observed a feature which is independent of 
rectification. No known line should be present at this 
wavelength, so we argue that it must be cosmetic. 

Special care has to be taken over the nebular con- 
tamination of the stellar spectra. The stars studied are 
located inside Hil regions, so the stellar spectra are 
contaminated by the nebular emission spectrum. It would 
be desirable to remove this nebular contribution, as it fills 
the cores of the Balmer H and Hei lines. However this is 
not easy even though we have long slit observations; the 
nebular emission has spatial structure that complicates its 
subtraction (especially for the Balmer H lines, the most 
important nebular contribution). If this is not correctly 
done an over or under subtraction will appear. After 
trying different possibilities we concluded that the best 
solution for this problem is not to subtract the emission 
lines and to ignore these regions in the final spectrum. 
For the Balmer H and Hei lines this is satisfactory, as 
emission lines are narrower than absorption lines. Nebular 
contributions could be more difficult to separate for metal 
lines, however the contamination of the stellar metal lines 
used for the abundance analysis due to nebular lines is 
negligible. 

The INT observations consist of the brightest three stars 
in the Trapezium cluster {9^ Ori A, C, D) together with 
the two nearby stars 9^ Ori A and B. Two standard stars 
were included in this set, 10 Lac and 15 Mon (09V and 
07V respectively). The other standard star, r Sco, was 
kindly provided by Dr. Gehren. This is a slow rotating 
B0.2V star perfect for a preliminary abundance analysis 
study. The spectrum was obtained with CASPEC, attached 
to the ESQ 3.6m telescope. The SNR of this spectrum 
is ^ 200-300 in the blue region and ^ 150 in the H^ region. 
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HD 


Name 


SpT 




Av 




Orion stars 


HD 37020 


0i Ori A 


B0.5V 


6.73 


1.89 


-3.4 


HD 37022 


Ori C 


07Vp 


5.12 


1.74 


-4.9 


HD 37023 


Ori D 


B0.5V 


6.71 


1.79 


-3.3 


HD 37041 


Ori A 


09V 


5.07 


1.12 


-4.3 


HD 37042 


Ori B 


B0.5V 


6.41 


0.73 


-2.6 


Reference stars 


HD 47839 


15 Mon 


07V 






-4.8 


HD 214680 


10 Lac 


09V 






-4.4 


HD 149438 


T ScO 


B0.2V 






-3.3 



Table 1. Identification, spectral type and photometric vi- 
sual data of the studied objects. Av and rriv values for Orion 
stars from Preibisch et al. II1999II . Mv values for these stars 
have been calculated considering a distance d ~ 450 ± 50 pc 
to the Orion nebula. Data for HD 214680 and HD 47839 are 
from Herrero et al. 1199211 . Photometric data for HD 149438 
are from Humphreys 11978t . Uncertainties in rriv, and Mv 
are 0.01, 0.03 and 0.3 respectively. 



HD 


4000 - 5000 


4600 - 5100 


H„ 


HD 37020 


170-280 


210 


220 


HD 37022 


170-330 


280 


290 


HD 37023 


220-450 


450 


320 


HD 37041 


200-275 


230 


210 


HD 37042 


160-220 


250 


230 


HD 214680 


140-190 


200 


160 


HD 47839 


89-120 


160 


260 



Table 2. SNR achieved for the different spectra for the three 
ranges observed with the int+ids. 



For the study of the spectral variability of Ori C, 
we have used feros spectra (some of them downloaded 
from the FEROS database and other kindly provided by 
Dr. Stahl). These observations were carried out with 
FEROS at the ESO 1.52m telescope in La Silla. The 
instrument is designed for high-dispersion spectroscopy 
with R - 48000 in the spectral range 3700 - 9200 A. The 
achieved SNR is 300 at about 4500 - 5000 A. Different 
phase observations have been used for the variability 
study (see Table [Tj). This study is presented in section [S| 

3. Determination of rotational velocities 

The analysis of stellar spectra makes use of a number 
of free parameters like the micro and macroturbulent 
velocities and the projected rotational velocity, usini. 
The last one has acquired particular importance in recent 
times because of the mixing that rotation may induce 
in the interior of massive stars (e.g. Maeder & Meynet, 
120001 Villamariz et al. . 1200211 . However, some methods to 
determine the rotational velocities do not distinguish be- 
tween rotation and other surface broadening mechanisms, 
like macroturbulence. 



Conventionally, v sin i values are based on linewidth 
measurements of individual features, mainly metal lines 
apparently free of blends. As the principal broaden- 
ing mechanism of these lines is stellar rotation, with 
sufficient resolution it is possible to determine v sin i 
from the FWHM of the line. Usually metal lines are 
used, however in cases of high rotational velocities or 
high temperatures metal lines appear blended or are 
very week. Therefore, in these cases, if the wsini is 
high, the whole He spectrum is used; if wsinz is not 
extremely high, only He i lines are used, as these lines 
are less affected by pressure broadening than Hell lines. 
However, the zjsini derived must be tested with some 
metal lines (if available), as we are not completely sure of 
rotation broadening dominating over pressure broadening. 

The Fourier method of determining v sin i is based 
on the fact that in the Fourier space convolutions 
transform into products and of the rotation, macrotur- 
bulence, natural and instrumental profiles (turbulence 
and instrumental assumed gaussian), only the rotation 
function has zeroes in its Fourier transform. These zeroes 
will appear in the total transform function, and Carroll 
H1933|l showed that the position of the zeroes are related 
to the usini. Actually the frequency of the first zero (cti) 
is related to the rotational velocity through: 

— V sin i CTi = 0.660 (1) 

The microturbulence also introduces zeros in the Fourier 
transform at high frequencies (Gray, I1973|l . This has 
to be taken into account for very low v sini (< 20 kms"'^). 

The main problems in the application of the Fourier 
method for early type stars are related to the quality 
of the observed spectra (i.e. spectral resolution and 
SNR). The lowest usini limit that can be determined is 
given by the spectral resolution (AA in A/pixel), as the 
sampling of the computational Fourier transform cannot 
be extended beyond 0.5/AA. The noise in the observed 
spectra transforms as white noise that obscure the first 
zero. 

The advantages of the Fourier method are that rotational 
broadening can be separated from other broadening 
mechanisms, and therefore metal lines as well as He lines 
can be used for the v sin i determination (even for low 
values of t)sini). This is very useful for fast rotating stars 
and spectral types earlier than 09, showing blended or 
very weak metal lines. 

We have tested the Fourier method with theoretical and 
observational cases, and it works well for massive hot 
stars. A paper with these results as well as determina- 
tions for a number of O star rotational velocities is in 
preparation (Simon-Di'az & Herrero, '2005|, figure |2I shows 
a typical example. For a recent application to A-stars see 
Royer et al. (2002,) . 
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Fig. 3. Four different synthetic lines generated with 
Fastwind have been convolved with a v sin i of 60 km s^^ 
and degraded to a SNR of 200. The different lines show 
the Fourier transform of the lines on a ctA baseline. The 
first zero is very close to the theoretical value for all the 
lines (even for the case of the He i line. 
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This 
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Other authors 




Fourier 


FWHM 


FWHM 


HD 37020 


55.0 ±0.6 


57 ±3 


56" 


HD 37022 


24 ±3 


40 ±11 


53'' 


HD 37023 


49.0 ±0.9 


51 ±3 


51° 


HD 37041 


131 ±4 


140 ± 11 


llO*" 


HD 37042 


34.0 ±0.5 


36 ±3 


10^ 


HD 214680 


30.0 ±0.8 


37 ±4 


35'' 


HD 47839 


67 ±4 


66 ±6 


67'' 


HD 149438 


< 13 


16 ±3 


5'', 19" 



Table 3. Projected rotational velocity derived from the 
Fourier and fwhm analyses. References from the literature are: 
" Simon-Di'az et al. I^HCT^ Howarth et al. [TM71 " McNamara 
& LarssonilSnS Schonberner et aLCnHHl " Killian et al. ilfljm 



The V sin i of our sample of stars has been deter- 
mined through the Fourier and fwhm methods. Results 
are presented in Tabled together with some values found 
in the literature, for comparison. The resolution in the IDS 
spectra is 0.23 A/pixel, so the lowest usini that could be 
detected through Fourier method is ^^20 kms~^. For 9^ 
Ori C FEROS spectra the resolution is ^0.03 A/pixel, so 
ijsini > 2 kms~^ for detection. For the r Sco spectrum, 
the resolution is 0.1 A/pixel, so the lowest detectable 
vsini is 8 kms^^. 

All the papers referenced in Table 13 use the fwhm 
method applied to the optical spectra of the stars except 
those by Howarth et al. I|1997|l . who use a cross corre- 
lation technique for lUE spectra, and Schonberner et al. 
l(TM5|l who compare the observed spectrum of r Sco with 



synthetic profiles. 

Comments on the individual stars' vsini determination, 
as well as the comparison between the values derived 
through Fourier and fwhm methods are presented in 
Section 14.11 Agreement between both methodologies is 
very good, however there are some interesting cases (see 
the study of 9^ Ori C in Sectional). 

4. Stellar parameters 

The analyses have been performed using the latest version 
of Fastwind (an acronym for Fast Analysis of STellar 
atmospheres with winds), a code which was originally 
described by Santolaya-Rey et al. H1997fl . See Puis et 
al. I|2UU5|I for the newest description of the code along 
with a discussion of comparisons with previous models 
and other spherical mass-losing codes. The latest version 
uses a more complete line-blanketing and a temperature 
correction method based on the energy balance of elec- 
trons. The technique used for the derivation of the stellar 
parameters is already standard and has been described 
elsewhere (Herrero et al.ESSl Rcpolust et al. O(M| . We 
only give here the main points. The analyses are based 
on visual fitting of hydrogen Balmer lines and Hel and 
Hell lines. Through the Hei/Heii ionization equilibrium, 
the effective temperature can be estimated; the wings 
of the Balmer lines are useful for the determination of 
the gravity and can give us some information about the 
stellar wind. 

The code also needs other parameters, such as micro- 
turbulence. He abundance and wind properties (mass 
loss, terminal velocity and /3 parameter). Actually 
wind properties are related through the Q parameter 
iQ = M /(vooR)'-^). 

Once the observed lines are fitted with the modeled 
ones, effective temperature, gravity. He abundance, mi- 
croturbulence and logQ are determined. The low density 
in the winds of the studied objects makes the spectrum 
insensitive to changes in Q, so that we can only determine 
upper values in most cases. Microturbulence has no effect 
on the H/He spectrum of early type stars of large gravity, 
as it has been shown by Villamariz & Herrero H2002|l . 
Therefore only effective temperatures, gravities and He 
abundances are determined for this step of the analysis. 
Of course, microturbulence is important for the derivation 
of metallic abundances and will be determined in the cor- 
responding section. The code also provides the emergent 
flux distribution and then mass, radius and luminosity 
can be calculated if is known (see Herrero et al. ll992|l . 

Errors in Tcff and log;; can be established generat- 
ing a microgrid around central values; visual comparisons 
between modelled lines and observations allow us to 
determine the range of possible values for these parame- 
ters. For further comments on the effects of varying the 
various physical parameters used in the analyses, their 
mutual interplay and their error bars see viz. Repolust 
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HD Name Teff log g R (Rq) M{Mq) log(L/L0) log Q 

± 1000 K ±0.1 dcx 



HD 37020 


Ori A 


30000 


4.0 


6.3 ± 0.9 


14 ± 5 


4.45 ± 0.13 


< 


-13.5 


HD 37022 


6^ Ori C 


39000 


4.1 


9.9 ± 1.4 


45 ± 16 


5.31 ± 0.13 






HD 37023 


Ori D 


32000 


4.2 


5.6 ± 0.8 


18 ± 6 


4.47 ± 0.13 


< 


-13.5 


HD 37041 


0^ Ori A 


35000 


4.1 


8.2 ± 1.1 


39 ± 14 


4.93 ± 0.13 


< 


-13.5 


HD 37042 


Ori B 


29000 


4.1 


4.5 ± 0.6 


9 ± 3 


4.11 ± 0.13 


< 


-13.5 


HD 214680 


10 Lac 


36000 


3.9 


8.3 ± 1.1 


20 ± 7 


5.02 ± 0.13 


< 


-13.5 


HD 47839 


15 Men 


40000 


4.1 


9.3 ± 1.3 


40 ± 14 


5.30 ± 0.13 




13.0 


HD 149438 


T Sco 


32000 


4.0 


5.6 ± 0.8 


11 ± 4 


4.47 ± 0.13 




13.0 



Table 4. Stellar parameters derived from Fastwind analysis. Only an upper limit for log Q can be derived for these stars. 
The microturbulences considered for the HHe analysis in each star are shown in the corresponding fitting plots. A normal value 
for the He abundance has been considered for all the stars (e = 0.09). 
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Fig. 4. Analysis of HD 37020 {6^ Ori A, B0.5V). Two models have been plotted for comparison with the observed 
spectrum. The adopted model has been represented with a solid line, while another model varying 1000 K in Toff has 
been plotted with a dashed line. Note the core of the stellar Balmer lines contaminated with emission from the nebula. 
Only the wings of the Balmer lines are used for the fitting. The narrow line that appears in the core of the He ii 4200 
line is a N ill absorption line; the red wing of the He 1 4922 line is contaminated by an O ii absorption line. The blue 
wing of the H^ line is affected by a cosmetic feature and is not used for the fitting 



et al. H2004|l . Herrero et al. 1)2002(1 . Villamariz & Herrero 
pir^ll . and Villamariz et al. plH^ . Errors in R, M and 
L are calculated considering the propagations of the 
uncertainties in Toff, log (7 and Mv. 

The fits of the synthetic Fastwind H and Hei-ii 



profiles to the observed ones are shown in Figures ^ to 
1101 The derived parameters for our sample of stars are 
shown in Table 0] corresponding to the best fits. 
Some comments on the individual analyses and the 
comparison between spectroscopic and evolutionary 
results are given in Sects . EtI and H31 
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Fig. 5. As Fig. 2|for HD 37023 {9^ Ori D, B0.5V). A variation of 0.1 dex in log 17 has been considered in this case. 



4.1. Comments on the Trapezium cluster stars 

HD 37020: This is an eclipsing binary discovered by 
Lohsen lfW75|l . The period is P = 65.43233 d, the 
magnitude range mv = 6.7-7.6 and the eclipse lasts 
~20 h. The primary is a B0.5V star which has many 
isolated metal lines. We have not found lines of the 
secondary star in our spectrum. This fact, together 
with the characteristics of the binary system described 
before, suggest that the secondary is a cooler, smaller 
and less luminous star, and that the change in the 
visual magnitude is due to the eclipse of the primary 
by its companion. Therefore the optical spectrum is 
mainly dominated by the B0.5V star. 

The Fourier method has been applied to some 11, 
Si iii-iv, N II and He I lines, deriving a v sin i = 55 ± 0.6 
kms~^. Metal and He I lines are in agreement. The 
V sin i derived through the linewidth measurement 
method is consistent with this value (see Table 

Figure ^ shows the good fit of the Fastwind 
profiles to the observed spectra for the parameters 
given in Table 01 (except for the forbidden component 
of He I A4471 which is not well reproduced throughout 
our analyses. Note also that the apparent bad fit of 



HeiiA4200 is due to the blend with the Niii line 
at the same wavelength). We are very close to the 
applicability limit of the He i-ii ionization equilibrium 
for deriving the Tes as HeiA4200 and Hei A4541 lines 
are very faint; however these lines are still sufficiently 
sensitive to changes in Toff and logg for deriving the 
stellar parameters accurately. The stellar parameters 
obtained here are in very good agreement with those 
obtained by Cunha & Lambert l|1992|l using the 
Stromgren index cq and the wings of H-,, with used 
Kurucz's fl979) LTE model atmospheres (this is not 
the case for the other two stars in common with these 
authors). A comparison of values is given in Table El 

HD 37022: This is the main ionising source of the Orion 
nebula. See section jsl for a detailed study of this star. 

HD 37023: This is the only star in the four Trapezium 
cluster stars {6^ Ori) without a binary companion 
(Preibisch et al. [TU^ . Robberto et al. ijMMjl find 
indications that this star is surrounded by a photoe- 
vaporated circumstellar disk. 

The Fourier method gives a w sin i = 49.0 ± 0.9 kms~^ 
for this star. Again, there is agreement with the 
linewidth measurement method. 
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Fig. m for HD 37041 {9^ Ori A, 09V). A variation of 1000 K in Toff along with 0.1 dex in log 5 has been 
in this case. Note that the broad wings of the Hei cannot be fitted correctly (see text) 



star Tetf (K) log g (dex) 

HD 37020 29970 / 30000 3.92 / 4.0 

HD 37023 32600 / 32000 4.70 / 4.2 

HD 37042 31600 / 29000 4.70 / 4.2 



Table 5. Comparison of stellar parameters for HD 37020, 
HD 37023, HD 37042. First values refer to Cunha & Lambert 
11992.) determinations, second values to this work. We see that 
there is excellent agreement for HD 37020, but poor agreement 
(specially for log g) for the other two stars. 



Figure El shows the fitting of the HHe lines. Observed 
He I lines are slightly broader than the theoretical 
ones. Table [S] shows the comparison between the 
stellar parameters we have derived and those by 
Cunha & Lambert H1002|l . In this case the agree- 
ment is not so good as for HD 37020, although the 
Teff are compatible, the log ^ they derived is very high. 

HD 37041: This is a single-lined spectroscopic binary 
(see Howarth et al. 1199711 . These authors find a single 
peak in the cross correlation function of the lUE 
spectrum, indicating that the spectrum of the primary 
is not seriously contaminated. 



Comparing the spectrum of this 09V star with 
that of the standard star 10 Lac (also classified as 
09V) we have found that there are no unblended 
metal lines due to its high rotational velocity (except 
SiivA4089, but it is in the blue wing of H^). A good 

V sin i determination has been possible using the 
Fourier method with the He i lines. A v sin i = 131 ± 4 
kms"""^ has been derived. We have used the Siiv 
line for checking the reliability of this value; a 

V sin i = 136 ± 5 km is obtained. He n does not give 
good results. A u sin i = 140 ± 11 kms^^ is derived 
using the fwhm method for the HeiA5015 line. 

Figure shows the fitting of the synthetic pro- 
files with the observed ones. A i;sini = 131 kms~^ 
has been considered for the H-He analysis. The wings 
of the He I lines cannot be well fitted; this might be 
explained by the presence of a companion. 

HD 37042: We have found good agreement between 
metal and Hei lines when using the Fourier method. 
The wsini derived is 34.0 ±0.5 kms~^. The value we 
obtain with the fwhm method is 36 ±3 kms"^, in 
agreement with the former one. 
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Fig. 7. As Fig. H for HD 37042 (9^ Ori B, B0.5V). For this star, the Hei lines fit better if a microturbulence of 10 



A very good fit of the observed and synthetic 
Fastwind profiles is obtained (see Figure [7|). Again 
the forbidden component of He i A4471 is too weak. 
For this star, the He i lines fit better if a microturbu- 
lence of 10 kms"^ is considered. 

This is the third star in common with Cunha & 
Lambert (see Table EJ; we also find for this star (as 
for HD 37023) that the stellar parameters derived by 
these authors are very different from ours (they obtain 
a very high log g and a Toft ^ 3000 K higher). 

4.2. Comments on the reference stars 

HD 21468O: This star has been previously considered 
as standard for stellar parameter determination 
through the H-He analysis (see Herrero et al. 120021 
and references therein). We have reanalysed here the 
spectrum of the star, obtained from a new observing 
run covering a larger spectral range (the HeiAA5015, 
5048, 6678 and HeiiA6683 lines can be used for the 
new analysis). We are using a new Fastwind version 
(Puis et al. 12005(1 ■ slightly different from that used 
by Herrero et al. The new analysis will be useful for 



checking the consistency of the latest version of the 
code in this stellar parameters regime. 

The projected rotational velocity of this star has been 
determined accurately by means of the fwhm method 
(wsini = 37±4 kms~^). The Fourier method applied 
to the iNT-fiDS spectrum gives a w sin i = 30 ± 0.8 
kms~^. This larger difference could be due to the fact 
that the v sin i is close to the computational Fourier 
transform limit (~20 kms~^ for this spectrum), or 
because the microturbulence is affecting the deter- 
mination of the rotational first zero in the Fourier 
transform fGrav ll975|l . 

The whole set of HHe lines are perfectly fitted 
with the Fastwind synthetic profiles (Figure 
The parameters derived by Herrero et al. (,2002|1 are 
rcff = 35500 K, log 3= 3.95 and e = 0.09. Our results 
are in agreement with those values. 

HD 149438: This star was selected for comparison in 
the stellar oxygen abundance analysis of the B0.5V 
Trapezium cluster stars (see section ^ . It has been 
studied elsewhere (Martin 0001 Przybilla & Butler 
I^HMl Kilian IBMl Schonberner et al. 119881 Becker 
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Fig. 8. As Fig. 2|for HD 214680 (10 Lac, 09V). A variation of 0.1 dex in log g has been considered in this case. 



& Butler iroggl Peters & Polidan [T^ Kane et al. 
ll98Ufl . Comparing isolated oxygen and silicon lines to 
a set of rotationally broadened profiles, Kilian et al. 
p991|l obtained a wsini = 19 kms~^ for this star. 
The Fourier method applied to the CASPEC spectrum 
only allow us to say that the vsini is lower than 13 



kmf 



this is because in this case the effect of the 



microturbulence over the broadening of the metal 
lines can be comparable with that produced by the 
rotation, so the first zero could be associated with the 
microturbulence instead of the vsini fGrav n"973|l . 



estimated a mass for the primary of M=34 M©, and 
suggested that the secondary has a spectral type 
09.5Vn (with very broad hues). 
We have used the spectrum of HD 214680 (09V) 
convolved with a high vsini (^350 kms~^), for 
recognizing lines in the spectrum of HD 47839 not 
contaminated by the secondary star; three metal lines 
were found. Using these lines (SiivAA4212, 4654 
and NiiiA 4379) a t;sini = 67±4 kms~^ has been 
determined. A similar value has been derived applying 
the FWHM method to the same lines (66 ±6 kms~^). 



A very good fit with the synthetic Fastwind 
profiles is obtained for this star (see Figure O. In 
this case the problem with the forbidden component 
of the He i A4471 can be clearly seen. Table [TTl (Sect. 
I6.4|l summarized the stellar parameters obtained in 
this and previous work. 

HD 47839: This star was selected for comparison with 
the main ionizing source in Orion (0^ Ori C). Some 
comments on the analysis of this star and the com- 
parison with 9^ Ori C are presented in section It 
was classified as 07V(f) by Walborn lfT!?7^ . Gies 
P993|l pointed for the first time that this star is a 
spectroscopic binary with a period of 25 years. He 



The fitting of the H and He lines for the stellar 
parameters show how the He I are contaminated by 
the secondary star lines. The spectroscopic derived 
mass (36 ± 9 Mq) is in very good agreement with 
the dynamical mass derived by Gies l)1993|l . Herrero 
ct al. iflM^ derived a Teg = 39500 K, log5=3.70 
and e = 0.07 for this star. Although we would ex- 
pect to obtain a lower Tgg due to the inclusion of 
line-blanketing effects, the value we have obtained is 
slightly higher because there is also a large difference 
between the iogg we derived (4.0) and the one by 
Herrero et al. (3.7). There is also another difference; 
we do not need a lower He abundance for fitting 
the He lines. This could be due to a binarity effect; 
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Fig. 9. As Fig.Elfor HD 149438 (r Sco, B0.2V). A variation of 0.5 dex in logQ has been considered in this case. The 
Ha line is not contaminated by nebular emission for this star and a more accurate value of log Q can be determined. 
The He i A 6678 is out of the observed range in the CASPEC spectrum of r Sco. 



when a composite spectrum is considered in a binary 
system, the lines can appear diluted or magnified due 
to the combination of the fluxes of the primary and 
the companion. If the system is out of eclipse the 
total flux will be higher than when the primary is 
observed isolated, so when the spectrum is normalized 
all the lines will appear diluted (and then a lower He 
abundance is needed to fit the He lines and a lower 
gravity is derived). 

4.3. Results of the stellar parameters study 

From the optical spectra of the Orion stars only upper 
limits for the Q parameter can be achieved. These 
estimations are based on the effect of the wind on the 
He II A4686 and H^ lines (the later one is contaminated 
by the nebular contribution). Some tests have shown 
that the other H and He lines are not affected by the 
uncertainties in logQ for the range of values considered, 
so the derived parameters will not be affected. 

Masses, radii and luminosities have been derived for 
all the studied targets (they are indicated in Table 01 
together with their uncertainties). The main source of un- 



certainty for these parameters is that associated with the 
absolute magnitude (except for very large uncertainties 
in logg). An error in Mv ~0.3 propagates to the mass, 
radius and logarithmic luminosity, giving uncertainties of 
~37%, 13% and 3% respectively. 

The stars have been plotted on the HR diagram in 
Figure Evolutionary tracks from Meynet & Maeder 
l)2()()3(l . corresponding to initial masses ranging from 9 to 
120 Mq and initial rotational velocities of kms~^ are 
also plotted. All stars are found in the Main Sequence 
close to the ZAMS, as it is expected because of their 
youth. Nevertheless, we can see the separation from 
the ZAMS increasing with luminosity, as pointed out 
by Herrero et al. H2()()4|l . The loci of the Orion stars 



is consistent with an isochrone at about 2.5±0.5 Myr, 
derived from the tracks with zero initial rotational 
velocity, which is slightly older than the upper limit given 
by Palla & Stabler (1 999. 2 Myr) and somewhat larger 
than other Trapezium age determinations found in the 
literature (e.g. Hillenbrand. 119971 <1 Myr). However, 
it has to be considered that, at large initial rotational 
velocities, the ZAMS is slightly shifted to the right and 
modifies the derived ages. Hence, until the role of the 
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Fig. 10. As Fig.Elfor HD 47839 (15 Mon, 07V). A variation of 1000 K in T^s has been considered in this case. 



initial rotational velocities is properly understood (for 
example, the distribution of initial rotational velocities in 
clusters), the use of isochones for massive stars in very 
young clusters should be regarded with special caution. 

A good agreement between gravities derived from 
the evolutionary tracks and those obtained from quan- 
titative analysis of the spectra is found (see Table EI). 
There is a trend for the most massive stars to have larger 
spectroscopic than evolutionary masses, but the number 
of objects is too small to draw any general conclusion. 



5. Modeling the main ionizing star of the Orion 
nebula 

HD 37022 [9^ Ori C, 07V) is the main ionizing source of 
the Orion nebula. We want to derive its stellar parameters 
as a first step to determine the effect of its ionizing flux on 
the photoionization of the surrounding nebula in a consis- 
tent way. Once these parameters are known, the spectral 
energy distribution can be modeled by using model atmo- 
sphere codes. In this way one of the inputs used in the 
photoionization codes will be fixed consistently. 



star 


A4poc (M0) 


A/ovol (Mq) 


log g spec 


log g ovoi 


0i Ori A 


14 


15 


4.0 


4.03 


Ori C 


45 


33 


4.1 


3.98 


Ori D 


18 


16 


4.2 


4.16 


Ori A 


39 


22 


4.1 


3.97 


Ori B 


9 


13 


4.1 


4.26 


10 Lac 


20 


25 


3.9 


4.02 


15 Mon 


40 


34 


4.1 


4.05 


r Sco 


11 


16 


4.0 


4.16 



Table 6. Comparison of masses and gravities derived from the 
evolutionary tracks and from the quantitative analysis of the 
spectra. The quoted log (?cvoi values are given corrected to two 
decimal places to be consistent with the corresponding evolu- 
tionary masses. Note however that these are not an indication 
of the precision of these values, which we consider to be 0.1 
dex. 



5.1. A historical review 

This star is known to have variable spectral features 
varying in phase or antiphase with a period of 15.422 
± 0.002 d (Stahl et al. IT^ Walborn & Nichols [TOMl 
Stahl et al. ll996|l . These variable features were discovered 
after Conti H1972|l showed for the first time that 6^ Ori 
C has a variable inverted P-Cygni profile in the Heii 
4686 line. Among them are Hq, emission, variability in 
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Fig. 11. HR diagram with the studied stars. Evolutionary 
tracks from Meynet & Maeder H2U(J3|I . Isocrones from 
Schaller et al. H1992|l . corresponding to 2, 2.5 and 3 Myr 

the equivalent width of some atmospheric and wind lines, 
and X-ray emission (Caillault et al. 119941 Gagne et al. 
I1997|l . Different explanations for this variability were 
postulated. The possibility of 9^ Ori C being a binary 
and this binarity explaining the spectral variability has 
been dismissed (Stahl et al. I1996|l . The variability has 
been associated with the rotation of the star. Stahl et 
al. (1996) proposed the presence of a dipolar magnetic 
field in 6^ Ori C, with the magnetic pole inclined 45° 
from the rotation axis (which is inclined 45° from the 
line of sight). The geometry of this system would imply 
alignment between magnetic pole and the line of sight 
at phase 0.5 and they would be perpendicular at phase 
0.0 (when maximum emission is found). Babel & 
Montmerle (|1997|l proposed the magnetically confined 
wind-shock model (MOWS) for explaining the variability 
in the star. According to this model, the radiatively 
line driven wind is confined by a dipolar magnetic field 
towards the magnetic equator of the system, generating 
a cold, dense disk due to the collision of material coming 
from both hemispheres, and a hot post-shock region. 

The wind characteristics of 9^ Ori C were deter- 
mined by Howard & Prinja lfT!W?|l and Stahl et al. ^T^M 
through UV spectrum studies. The former derived a mass 
loss rate of 4x10""^ Mq yr"-'^, the latter through the 
absorption in CiV lines, determined a terminal velocity 
somewhat greater than 2500 kms~^. It is common that 
07V stars have stellar winds; what is not so common is 
the detection of the presence of magnetic fields in O stars. 
Donati et al. H2002|l succeeded in the detection of Zeeman 
features in the spectrum of 9^ Ori C through spec- 
tropolarimetric observations with the Anglo- Australian 
Telescope. They detected variability in the Stokes V 
profiles of some photosphcric metal lines. This variability 
has a coherent modulation with the period derived from 
other variable features. However, the geometry derived 



was in contradiction with that from Stahl ct al. (CHMI), 
with the magnetic pole aligned with the line of sight 
at phase (they found a maximum in the longitudinal 
component of the magnetic field at this phase) . 

5.2. Preliminary study of the spectrum 

Having spectral variability, it is very important to 
understand the cause of this variability and to deter- 
mine which lines are reliable for the stellar atmosphere 
modeling before comparing synthetic and observed 
H - He profiles. Preliminary work with the int-|-ids 
spectrum showed that a better spectral resolution was 
needed to apply the Fourier method for obtaining the 
vsini. This spectrum did not allow us to have neither 
a good enough sampling (AA), nor to carry on a study 
of the variability, so we decided to use some feros 
spectra with better quality and covering all variability 
phases (see TableO, that are available in the ESO archive. 

Through the study of the feros spectra observed 
in the different phases we have found all the variable 
spectral features described in Stahl et al. H1996|l and some 
new ones: 

— He II 4686: Broad emission appears in the blue wing 
of the line (the so-called inverted P-Cygni profile, with 
maximum at ~ 0) . Broad emission is also present 
in the red wing (maximum at (/)^0.5, minimum at 
0~O.8). 

— Balmer lines: These lines are affected by the same 
broad emission features than those in Hen A4686. The 
emission features are stronger in Hq, and H^ than in 
the other Balmer lines 

— Metal and He i-ii lines: All the line strength varies 
in phase, being larger at ~ 0. 



Name 


Date 


A/,7r>-2.400.000.5 





f07341+51 


16/10/98 


51102.31 


0.180 


f85221 


26/07/99 


51385.43 


0.539 


f03551+61 


08/10/98 


51094.28 


0.659 


f04711+21 


10/10/98 


51096.39 


0.796 


fl5241 


28/11/98 


51145.37 


0.972 



Table 7. feros spectra used for the study of the spectral 
variability of 9^ Ori C. All spectra have been downloaded from 
the ESO-FEROS database except 185221, kindly provided by O. 
Stahl. The diflferent phases have been calculated from M,JDo - 
2.400.000,5 = 48832,5 (Stahl et aL lTOga and P= 15.422 days 



This variability can be easily explained considering the 
model proposed by Stahl et al. (1996) and developed by 
Babel & Montmerle (1097). According to this model we 
would have an 07V star with a disc. The disc is produced 
by the confinement of the wind by a dipolar magnetic 
field through the magnetic equator. We will have a cool 
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Fig. 12. The most representative phases have been selected for showing the variability of the Hell A4686, and Hq, 
lines in the 9^ Ori C spectrum. The variable feature associated with the inverted P-Cygni profile discovered by Conti 
in the HeiiA4686 line is also present in the hydrogen Balmer lines. Another emission feature can be clearly seen in 
the red wing of these lines at phase c^O.S. The narrow emission features in the Balmer lines arc nebular lines. 



disc with material falling back to the stellar surface. If 
we consider that at phase the cool disc is seen edge-on, 
the blue shifted emission appearing in He li A4686, H^ and 
the other hydrogen Balmer lines can be explained as stel- 
lar photons absorbed and reemited with a doppler shift 
corresponding to the velocity of the disc material falling 
to the surface of the star (in a process similar to what 
occurs in a stellar wind but with blue-shifted emission 
and red-shifted absorption). As density in the disc is very 
high then a strong blue-shifted emission will appear. At 
phase 0.5, when the disc is seen face-on the blue-shifted 
emission disappears. The emission appearing in the red 
wings of the former lines could be explained as the effect 
of the scattering of stellar photons by the wind material 
confined by the magnetic field and accreting to the disc 
(see figure lT^ . The disc will also have continuum emission 
that will affect the total continuum flux received from the 
star. This effect will be maximum when the disc is face-on 
because the emitting region is larger at this phase. The 
variability observed in Hei, Heii and metal lines (except 
for the emission in HeiA4686) is only a consequence of 
this effect. As the total continuum flux is varying with 
the phase, the normalized spectrum will be affected. All 
absorption lines will be artificially weaker when the disc 
continuum emission is at maximum. This variability in the 
lines can allow us to estimate the amount of continuum 
flux that comes from the disc, and then the visual magni- 
tude variability. By assuming that at phase 0.0 the lines 
are not affected by the continuum emission from the disc, 
we have scaled the spectra at the other phases to fit in 
the former spectrum; the scaling factor will be related to 
the ratio of visual fluxes (i.e. the difference in magnitudes) 
between the stellar component and the stellar -|-disc contri- 
bution. This can be seen in Figure^] From this study we 
would expect a variability of ~ 0.1 6 magnitudes. Kukarkin 
et al. I)1981|l . in their catalogue of suspected variable stars. 
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Fig. 13. Magnitude variability expected from the study of 
the He ii A5411 line. The other He i-ii lines follow a similar 
behavior. The presence of a disk could be responsible of 
this variability (see text for details). The solid line corre- 
sponds to a sinusoidal curve with a maximum change in 
mv of 0.16, presented for comparison. 

found a change in ra^ of 0.06 magnitudes (5.10 - 5.16) for 
9"^ Ori C. Hipparcos has also classified this star as variable; 
although Hipparcos data do not show a clear pattern, the 
median magnitude in Hipparcos system is 4.61 mag and 
the variability of this magnitude varies between 4.56 and 
4.70, that is agreement with our study. 

5.3. Determination of the v sin i of 9^ Ori C 

It is very important to have a good determination of the 
projected rotational velocity of this star as it is supposed 
that the spectral variability of 9^ Ori C is related with 
its rotation. The derived usini should be independent 
on the phase and should be coherent with a period of 
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~15.4 days. It is shown in section [^| that the profiles of 
metal as well as He i and He ii lines are dependent on the 
phase, however this is only an artificial dependence due 
to the presence of the disk. Once the spectra of different 
phases are corrected from the effect of dilution by the 
disc continuum emission, all the metal, He i and He ii 
lines are independent on the phase (except those related 
with the disc, see Figure IT^ . 

The Fourier method allows us to separate pure ro- 
tational broadening from other broadening mechanisms 
affecting the shape of the lines. We have used this method 
with some metal lines at phase ~ (see tabic IHl . A 
usini = 24±3 kms~^ has been derived. Figure ITU shows 
an example of the application of the Fourier method in 
the determination of the vs'mi of 6^ Ori C. 

Table 7 offers a comparison of v sin i values obtained with 
the Fourier and fwhm methods. We see that the Fourier 
method gives more consistent values for all considered 
lines. Differences within the FWHM method may reach a 
factor of 2 and, in fact, we had problems when trying to 
fit the profile of some of the lines with a gaussian profile 
for measuring their linewidth. 

The derived value for v sin i through the Fourier method 
is also more coherent with a 07V star rotating with a 
15.422 days period. Considering R^U Rq, the upper 
limit for vsini is ~35 kms~^and therefore the inclination 
of the rotational pole is i ~ 45, in agreement with previous 
results (see next section). 



Line 


V sin i 


(kms^^) 




Fourier 


FWHM 


C in 4056 


23 


42 ± 10 


Niii 3998 


23 


28 ± 7 


Niv 4057 


20 


33 ± 7 


Niv 5200 


20 


32 ± 10 


Si IV 4089 


23 


62 ± 14 


Si IV 4654 


30 


47 ± 10 


Oiii 4081 


21 


52 ± 10 


O III 4376 


26 


31 ± 4 


O III 4435 


25 


32 ± 5 


O III 4454 


27 


34 ± 4 


Oiii 5592 


25 


45 ± 6 



Table 8. Projected rotational velocities derived through 
Fourier and fwhm methods for some metal lines present in 
the spectrum of 9^ Ori C. Values were derived at phase 0.972 
(see explanation in text). 



5.4. Modeling of 6^ Ori C 

Once the spectral variability of 6^ Ori C is understood 
and its vsini has been derived, we can proceed to model 
its stellar atmosphere and wind through the observed 
spectrum of the star. The lines used for this analysis are 
shown in Figure 1151 basically they are the ones used in 




0.1 1.0 10.0 

a (cicle/A) 

Fig. 14. Fourier analysis of the N iv A4057 line in 6^ Ori 
C. 



the other analyses plus Hei 5875 A. Some of the lines 
are contaminated (see Sect. 15. 2() . so this must be taken 
into account. The H^ and H^ lines have been selected 
as the most reliable lines for deriving logg (they are less 
contaminated than Hp and Hq); The whole set of Hei-il 
lines has been considered except HeiiA4686; however 
it must be taken into account that the strength of all 
these lines vary with the phase. Phase 0.972 will be used 
for determining the stellar parameters (the effect of the 
continuum emission of the disc is smaller at this phase, 
see Sect. l5.2|l . 

Although our study has shown that the rotational 
velocity (derived from Fourier analysis) is ~ 24 kms~^, 
when this broadening is considered all synthetic lines 
appear narrower than the observed ones. We have tried 
to solve the problem by means of a different vsini value, 
however it does not work because then the shape of the 
synthetic profiles do not fit with the observed ones (the 
cores of the modeled profiles are too wide when the FWHM 
of the He lines is fitted). An extra-broadening mechanism 
has to be included. When a gaussian macroturbulent 
broadening (Grav. 11^73;) is used, the fit clearly improves, 
however in this case the Hi and Hen lines cannot be 
fitted simultaneously with the He I ones; when the former 
are fitted (for a T'eff = 39000 K), some of the synthetic 
lines in the latter appear stronger and narrower than 
observed. A better fitting for the Hei lines is obtained if 
a higher Tcs if considered, but then the synthetic Hen 
lines appear too strong. There is no way to solve this 
problem in this region of the parameter space either by 
changing the rotational velocity, the macroturbulence or 
the microturbulence. 

In section 14.21 figure EH shows that the fitting of 
the He i-n lines could follow a similar behavior in the case 
of the 07V star HD 47839 (selected as reference star). A 
variation of 1000 K in the effective temperature change 
strongly the strength of the Hex lines. For this spectrum 
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4540 
X (A) 



4685 
X (A) 



4695 5402 



5412 
X (A) 



5422 6670 



6682 
X (A) 



6695 



Fig. 15. HHe analysis of HD 37022 (9^ On C, 07Vp) at phase (/) = 0.972. The sohd black line corresponds to a model 
with Teff, = 39000 K, logg=:4.1 e = 0.09, ^t = 5 kms^^ and logQ = -13.0; the dashed blue line corresponds to a model 
with Toff =40000 K (same remaining parameters). A wsini = 24 kms~^ and a macroturbulent velocity of 60 kms^^ 
have been considered. The observed lines strength has been enlarged by a factor 1.1 to correct for the effect of the 
disc continuum emission (see text for explanation). 



also occurs that when the He I lines appear fitted, the 
Hen lines are slightly stronger than observed, and if 
the Hell lines are fitted, the He I lines are stronger than 
observed. 

Puis et al. H2005|l have shown that there is a dis- 
crepancy for the He I singlets between the synthetic 
Fastwind and CMFGEN lines for effective temperatures 
between 36000 and 41000 K (being the CMFGEN profiles 
shallower). Therefore the Hei triplet system can be 
considered more reliable (i.e. the HeiA4471 line). 
Knowing this discrepancy, we have considered the 
HelA4471 for the fitting with Fastwind synthetic pro- 
files. Our best model corresponds to Tes = 39000 ± 1000 
K and log g= 4.1 dex. From these values, the model flux 
and with Mv=-4.9±0.3, a stellar radius 7?= 10.6 ±1.5 
Rq is derived, and then an inclination of the rotational 
axis of 1 = 44 ±12°. This value is in agreement with 
previous independent studies (Stahl et al. 119961 Donati 
et al. I20()2|l , although our derived v sin i is more reliable 
and our radius is not obtained from a SpT - R calibration 
but is a result of the spectral analysis of the star. 



6. Oxygen abundances in the Trapezium cluster 
B0.5V stars 

Three of the stars studied in Orion are perfect targets for 
a stellar abundance analysis as they have many narrow 
unblended fines. These are HD 37020, HD 37023 and 
HD 37042, three B0.5V stars. Cunha & Lambert 
T994I presented carbon, nitrogen, oxygen and silicon 
abundances from LTE and NLTE analyses for these 
stars. For comparison purposes, a similar analysis has 
been done for r Sco, a B0.2V star (Walborn & Fitzpatrick, 
I1990|) with very low v sin i. Stellar abundances for this star 
have been derived elsewhere in the literature (Hardorp 
& Scholz lWn Kane et al. 1980; Peters & PofidanHnHl 
Schonberner et al. I1988i Becker & Butler ,.1988. Kilian et 
al. lTMTl Martin 12001 see Table [TT|l . The other two Orion 
stars have been ruled out: HD 37041 has a relatively 
high projected rotational velocity, so metallic lines are 
broadened and then appear blended; HD 37022, being 
an 07V star, does not have enough oxygen lines for an 
accurate abundance analysis. 
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We have therefore derived oxygen abundances for 
the Trapezium cluster B0.5V stars for comparing them 
with the M42 nebular abundances obtained by Esteban 
et al. (HXMJ. 

6.1. Line identification and measurement of equivalent 
widths 

We make use of the classical method of curve of growth to 
determine the oxygen abundances. When this methodol- 
ogy is used it is important to remove all lines that appear 
blended. The spectrum of r Sco (a star with similar spec- 
tral type than our targets in Orion and a very low vsini) 
has been used for an identification of the O ii lines present 
in the spectra of the Orion stars. The whole set of lines is 
shown in Table [T^ divided into multiplets, together with 
them log gf values (basically they are taken from the NIST 
database). 

The equivalent widths of all the O ii lines listed in Table 
1131 have been measured for the four stars, however only 
our set of suitable lines (see below) is shown in Table 
To measure the equivalent widths we use our own software 
developed in IDL. A least squares profile fitting procedure 
was used, with Gaussian profiles fitting the line and poly- 
nomials of degree one or two to fit the local continuum. 
Errors in the measurements due to the uncertainty in the 
position of the local continuum (estimated as ±1/SNR, 
Villamariz et al. I2()()2|l have also been considered. The es- 
timated value of the uncertainty in the measurement of 
the EWa is ~5 mA and ~10 mA for some problematic 
lines. 

6.2. Line selection 

Some of the lines that appear unblended in the spectrum 
of T Sco cannot be used in the analyses of the other stars; 
they have larger rotational broadening and then appear 
blended (or lie in the wings of H or He lines) . Once a first 
set of unblended lines was selected, a preliminary abun- 
dance analysis was done separately for each multiplet. 
In this way the dispersion in the line abundances for the 
zero slope value of the microturbulence are minimized as 
all the lines in a multiplet are formed in the same region 
in the stellar atmosphere. Problematic lines, errors in the 
measurement of the equivalent widths or artificial trends 
can be then identified; such lines will be removed in the 
global analysis (e.g. this is the case for the OiiA4414 
line, an isolated line that gives systematically lower 
abundances) . 

The set of suitable lines finally used in the abun- 
dance analyses is presented in Table |51 They are 
lines coming from transitions between configurations 
2p2 (3p)3s-2p2 (3p)3p (NIST multiplets 64, 65 and 
72) and 2p2 (3p) 3p- 2p2 (3p) 3d (NIST multiplets 90, 
148 and 130). We have ruled out the lines that do not 
follow the general trend. We have found that lines from 



multiplets 99, 118, 161, 188, 172 give systematically lower 
abundances. This effect can be due to the definition of the 
oxygen model atom we have used, or can be associated 
with the log gf values. Some preliminary comparisons 
with TLUSTY analyses have shown that the difference in 
the line abundances is also present for the case of r Sco. 

6.3. Chemical analysis 

For each star we proceed as follows. A grid of 16 
Fastwind models combining four abundances and four 
values of microturbulence is calculated. In this way, 
the curves of growth for each line can be constructed 
by plotting the theoretical equivalent width for each 
abundance and microturbulence versus the abundance 
(see Figure El for the case of the On A 4414 line in 
HD 37042). 

Through the observed equivalent width and its error one 
abundance (and its derived uncertainty) can be derived 
for each line and microturbulence. The individual line 
abundances are dependent on the microturbulence. The 
microturbulence value that minimizes the dependence of 
the line abundances on the line strength in the loge- EW 
diagrams will be the microturbulence we are looking for. 
These diagrams can be also used as a diagnostic tool to 
check the reliability of the different lines for the abun- 
dance determination (see Section f6.2|l . Figure [T7I shows 
the log e- EW diagrams for two different microturbulences 
in the study of HD 37023. 




I I I I 

8,35 8,50 8,65 8,80 



log (0/H) +12 

Fig. 16. Example of curve of grow for the line OiiA4414 
in the star HD 37042. A grid of 16 models has been con- 
sidered (4 microturbulences and 4 oxygen abundances). 
The observed EW of the line and its uncertainty are plot- 
ted as horizontal lines. Two examples of abundances (and 
their uncertainties) derived for microturbulences 7 and 9 
kms~^ are plotted as vertical lines. 

The microturbulence derived from the zero slope for each 
star is presented in Table ITUl Uncertainties in the micro- 
turbulence are obtained considering the errors derived 
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Line 


EWa 


e{Or 


Ae(0) 


EWa 


£(0) 




EWa 


e(0) 


Ado) 


EWa 


£(0) 


Ae(0) 




HD 149438 (^t 


= 8.7) 


HD 37020 


t = 6.5) 


HD 37023 (5 


t = 7.4) 


HD 37042 ({ 


t = 5.5) 


Toff = 


32000 K 


log 3 = 4.0 


Toff = 


30000 K 


log g = 4.0 


Toff = 


32000 K 


log 3 = 4.2 


Toff = 


29000 K 


log g = 4.2 


On 4638 


85 


8.75 


0.04 


— 


— 


— 


— 


— 


— 


— 


— 


— 


On 4641 


127 


8.65 


0.04 






— 






— 






— 


On 4661 


90 


8.70 


0.04 


80 




— 


76 


8.60 


0.06 


88 


8.56 


0.06 


On 4676 


79 


8.71 


0.05 


82 


8.59 


0.06 


65 


8.59 


0.06 


82 


8.63 


0.06 


On 4317 


92 


8.69 


0.03 


107 


8.67 


0.04 


78 


8.59 


0.05 


98 


8.60 


0.07 


O TT A'il Q 


86 


8.71 


04 


100 


8.69 


05 


75 


8.63 


05 


95 


8.66 


07 


On 4366 


80 


8.61 


0.04 


96 


8.62 


0.05 


77 


8.62 


0.05 


103 


8.62 


0.08 


On 4416 


100 


8.62 


0.04 


120 


8.65 


0.05 


93 


8.62 


0.05 


117 


8.70 


0.11 


On 4452 


37 


8.67 


0.07 


48 


8.64 


0.07 


32 


8.58 


0.09 


60 


8.67 


0.10 


On 4072 








118 


8.61 


0.12 


99 


8.57 


0.09 


119 


8.71 


0.10 


On 4076 








133 


8.59 


0.12 


114 


8.57 


0.09 








On 4078 








46 


8.67 


0.17 


33 


8.59 


0.09 


48 


8.68 


0.10 


On 4941 


38 


8.64 


0.06 


50 


8.68 


0.07 


35 


8.60 


0.09 


43 


8.60 


0.08 


On 4943 


60 


8.62 


0.04 


67 


8.58 


0.06 


55 


8.59 


0.07 


56 


8.48 


0.07 


On 4956 


17 


8.73 


0.12 


20 


8.67 


0.14 


13 


8.58 


0.20 


22 


8.71 


0.13 


On 4891 


24 


8.75 


0.10 


27 


8.68 


0.11 


18 


8.59 


0.15 


27 


8.65 


0.12 


On 4906 


34 


8.64 


0.07 


42 


8.63 


0.08 


31 


8.59 


0.10 


47 


8.72 


0.08 



£(0) = 8.70 ± 0.10 e(0) = 8.65 ± 0.10 £(0) = 8.59 ± 0.10 e(0) = 8.64 ± 0.10 



" €(0) = log(0/H) + 12 

Table 9. Equivalent widths and derived line abundances for the set of On lines used in our analysis. Line abundances refer 
to the microturbulence given in brackets for each star (^t in kms~^). Uncertainties in the line abundances come from the 
propagation of the uncertainties of the equivalent width measurements (see text). Some of the O ll lines of the Orion stars have 
not been used as they appear blended. O ll AA 4072, 4076 and 4078 lines were ruled out in the analysis of r Sco due to the poor 
quality of the CASPEC spectrum in this region. Final oxygen abundances for each star have been calculated through a weighted 
mean of the linear values. Errors represent the statistical deviation for these mean values. 
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Fig. 17. Example of loge- EW diagrams in the study of 
HD 37023. 

for the slope in a linear fit of the data (due to errors in 
the individual line abundances). This step also allows us 
to estimate the contribution of the uncertainty in the 
microturbulence to the the total oxygen abundance. This 
uncertainty depends mainly on the quality of the spectra 
(in this case is ~ 0.06 dex). 

Abundance values for each line as well as their un- 
certainties are obtained for that microturbulence (see 
Table El). The final abundance value is calculated through 



a weighted mean of the linear individual line abundances 
(lO*^'"^^), and its uncertainty (ctw) is that associated with 
this mean. Final values for the total oxygen abundances 
for each star are shown in Table EH The final uncertainty 
in the oxygen abundance takes into account four different 
sources of errors: those associated with the statistical 
analysis, those derived from the error in the determined 
microturbulence and finally those refered to the uncer- 
tainty in the stellar parameters and to atomic data. All 
these sources of error are added quadratically for deriving 
the final abundance uncertainty (see Villamariz et al. 

12001. 



6.4. Oxygen abundance in the stardard star r Sco 

Table sumarizes the oxygen abundances appeared in 
the literature for t Sco. Our derived value is compatible 
with previous results but a little higher than most of 
them. This difference can be easily explained taking 
into account that for this range of stellar parameters, 
the oxygen abundance derived through O ii lines is very 
sensitive to a change in T^g and log 5 (the lines become 
fainter when a higher Tcff is considered and then the 
derived oxygen abundance is higher). The uncertainty in 
the oxygen abundance due to a change of ~ 1000 K in Tcs 
is ~0.08 dex. This effect is considered in the uncertainty 
for the given values, however the central value will be 
slightly dependent on the derived stellar parameters. We 
have considered two Fastwind models with different Tes 
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Object 


N 


?t (km/s) 


log(0/H) + 12 


This work 


HD 149438 


14 


8.1 ± 1.2 


8.70 ± 0.10 


HD 37020 


13 


6.8 ± 1.3 


8.65 ± 0.10 


HD 37023 


15 


6.3 ± 3.0 


8.59 ± 0.10 


HD 37042 


13 


4.9 ± 2.0 


8.64 ± 0.10 






Cunha & 


Lambert (LTE) 


HD 37020 


6 


7.0 


8.92 ± 0.05 


HD 37023 


7 


7.5 


8.76 ± 0.04 


HD 37042 


7 


6.0 


8.97 ± 0.04 






Cunha & 


Lambert (NLTE) 


HD 37020 


7 


5.0 


8.83 ± 0.12 


HD 37023 


6 


7.0 


8.87 ± 0.08 


HD 37042 


6 


6.0 


8.85 ± 0.06 


Estcban ct al. 


Gas 






8.65 ± 0.03 


Gas + dust 






8.73 ± 0.03 



Table 10. Oxygen abundances for the three B0.5V stars in- 
side Orion Nebula and the reference star r Sco. Oxygen LTE 
and NLTE abundances derived by Cunha & Lambert (I1994ll 
for the Orion stars as well as those calculated by Esteban et 
al. lUTOll for the nebula are also presented for comparison. 



(32000, 32500 K) and the same log 5 obtaining oxygen 
abundances of 8.70 and 8.74 dex respectively. 

The derived abundance is also very dependent on 
the microturbulence considered (specially if lines with 
high equivalent width are used). We have taken into 
account these dependencies in our uncertainties. 



Work 


Toff (K) 


logs 


e(0) 


?t (kms-i) 


Hardorp & Scholz 


30900 


4.05 


8.70 




Kane et al. 


30300 


3.95 


8.6 




Peters & Polidan 


31500 


4.3 


8.65 ±0.26 


5 


Schonberner et al. 


33000 


4.15 


8.5 


5 


Becker & Butler 


33000 


4.15 


8.70 ±0.19 


5 








8.58±0.19 


10 


Kilian et al. 


31400 


4.24 


8.61 ±0.12 


3 


Martin 


30000 


3.75 


8.58±0.17 


7 


This work 


32000 


4.0 


8.70 ±0.10 


8 



Table 11. Comparison of stellar parameters and abundances 
derived for r Sco in previous studies found in the literature 
and in this work. 



6.5. Comparing with nebular and previous stellar 
results 

The oxygen abundances derived for the Orion stars are 
compatible within the errors (see Table I10() . HD 37023 
has a slightly lower abundance (but still compatible with 
the other abundances). In section 1^1 we have seen that 
the fitting of the H and He lines is not so good as for 
the other B0.5V stars (observed lines appear slightly 
broader). 



The oxygen abundances in the Orion stars are systemati- 
cally lower than that derived for r Sco. 

Esteban et al. H2004|l have recently published a reap- 
praisal of the chemical composition of the Orion nebula. 
They derived a total oxygen gas-phase abundance 
e(0) = 8.65 ± 0.03. However, some oxygen is expected 
to be depleted onto dust grains in ionized nebula, so 
the total gas-|-dust oxygen abundance should take this 
depletion into account. In a previous work (Esteban et 
al. I1998|l . these authors estimate the depletion onto dust 
grains by comparing Si and Fe nebular abundances with 
those obtained by Cunha & Lambert fT?)M| for B stars 
in the Orion association, assuming a certain composition 
for the main dust molecules. Taking into account this 
correction, the final gas-|-dust oxygen abundance Esteban 
et al. pIMjl propose is e(0) = 8.73 ± 0.03 (an oxygen 
abundance correction for dust ~ 0.08 is applied). 

Our stellar results are compatible with those ob- 
tained by Esteban et al. H2004|l for the gas phase, however 
the dust-l-gas corrected abundance seems to be too 
high compared with our derived stellar abundances, al- 
though still marginally consistent inside the uncertainties. 

Our oxygen abundances are systematically lower than 
the NLTE abundances by Cunha & Lambert H1994|l . In 
that paper the former authors comment that the LTE 
abundances are slightly more reliable than the NLTE 
abundances they present. The difference between their 
LTE abundances and our results are even higher. These 
differences can be in part associated with the differences 
in the derived stellar parameters for these stars. The Tgg 
and log g obtained by these authors are higher than ours 
(see Table El, so the oxygen abundances they derive are 
obviously higher. 

This difference in the effective temperatures may also 
affect the derived silicon abundances (used by Esteban et 
al. 119981 for estimating the oxygen depletion). Cunha & 
Lambert H1994(l . derive their silicon abundances by using 
3 Sim lines. Preliminary silicon analysis by our group 
has shown that a difference of ^ 1000 K in Tgg can shift 
the Sim abundances up to 0.2 dex (deriving a higher 
abundance for the lower Toff). 

Alternatively, our result could suggest that the molecules 
that Esteban et al. used to estimate the O dust depletion 
in Orion cannot be present in large amounts in this 
nebula. Consequently, Si, Mg and Fe (the refractory 
elements being the main constituents of those molecules) 
have to form other molecules without oxygen. 



7. Conclusions 

By means of a detailed spectroscopic analysis of the 
optical spectra of the Trapezium cluster stars, stellar 
parameters and oxygen abundances have been derived. 
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Label A (A) log gf Notes 



2s^ 2p^ (3P) 3p - 28^ 2p^ (3P) 3d 



NIST 90 (4Do-4F) 


O II 4069 


4069.623 


0.149 + Cm 




4069.882 


0.344 


On 4072 


4072.153 


0.552 


O II 4076 


4075.862 


0.693 +Cii 


O II 4078 


4078.842 


-0.284 


NIST 118 (2Do-2F) 


On 4699 


4699.218 


0.269 




(4699.011 


0.418) NIST 172 




(4698.437 


-0.883) NIST 172 


On 4705 


4705.346 


0.477 


On 4741 


4741.704 


-0.987 Very weak 


NIST 130 (4So-4P) 


On 4891 


4890.856 


-0.436 


O n 4906 


4906.830 


-0.160 


NIST 148 (2Po-2D) 


On 4941 


4941.072 


-0.054 


On 4943 


4943.005 


0.239 


On 4956 


4955.707 


-0.573 Very weak 


2s^ 2p^ (3P) 3s - 2s^ 


2p^ (3P) 3p 


NIST 64 (4P-4Do) 


On 4638 


4638.856 


-0.332 -1- Cn H- Sim 


On 4641 


4641.810 


0.054 -1- Nni 


On 4650 


4649.135 


0.308 


On 4651 


4650.638 


-0.362 


On 4661 


4661.632 


-0.278 


On 4673 


4673.733 


-1.088 


On 4676 


4676.235 


-0.394 


On 4696 


4696.352 


-1.380 


NIST 65 (4P-4Po) 


On 4317 


4317.139 


-0.386 




4317.696 




On 4319 


4319.630 


-0.380 -l-Nni 




4319.866 


-0.502 


On 4366 


4366.895 


-0.348 




4366.530 




NIST 72 (2P-2Do) 


On 4414 


4414.899 


0.172 




4414.456 


-1.483 


On 4416 


4416.975 


-0.076 -1- Si IV 


On 4452 


4452.378 


-0.789 



Table 12. Preliminary set of O ll lines selected for the analy- 
sis, divided by multiplets. The spectrum of the low v sin i star 
T Sco, has been used to identify the lines, log values are from 
the NIST database. 



Projected rotational velocities has been obtained through 
Fourier method. This method has been extensively used 
for late type stars, but not widely applied to early type 
stars. Our results show this method to be very useful for 
distinguishing between rotational broadening and another 
broadening mechanisms that can be present in early type 
stars (i.e. macroturbulence). The agreement is very good 
when comparing with results from the line- width method. 
The Fourier method applied to the high resolution 9^ 
Ori C FEROS spectra, allow us to derive a very accurate 
V sin i that is in agreement with the period of variability 
of some spectral features in 9^ Ori C. 



Label 


A (A) 


log gf Notes 


2s^ 2p^ 1 


;iD) 3s - 2s^ 


2p^ (ID) 3p 


NIST 99 (2D-2Fo) 


On 4590 
On 4596 


4590.974 
4595.957 
4596.177 


0.350 
-1.032 
0.200 


2s^ 2p^ 1 


;iD) 3p - 2s^ 


2p^ (ID) 3d 


NIST 161 (2Fo-2G) 


On 4185 
On 4189 


4185.440 
4189.788 
4189.581 


0.604 -l-Ciii 
0.716 
-0.828 


NIST 188 (2Po-2P) 


On 4691 
On 4701 


4690.888 
4691.419 
4701.179 
4701.712 
4700.441 


-0.610 
-0.309 
0.088 
-0.611 
-3.298 


NIST 172 (2Do-2F) 


O II 4703 


4703.161 


0.263 



Table 13. (Cont) Preliminary set of On lines selected for the 
analysis, divided by multiplets. The spectrum of the low v sini 
star r Sco, has been used to identify the lines, log g/ values are 
from the NIST database. 



Stellar parameters and their uncertainties have been 
derived for the studied stars using H , He i and He Ii lines 
and Fastwind code. 

The presence of many O Ii lines in the optical spec- 
trum of three B0.5V Orion stars has allowed us to work 
on a very detailed abundance analysis using the curve 
of growth method. This analysis has been performed 
through a careful selection of suitable lines from a 
previous study of the different O ii multiplets. In this 
way, the dispersion in the line abundances is reduced, 
and the final abundance value derived is very precise. 

The derived oxygen abundances in the Orion stars 
are in agreement with the nebular gas-phase abundances 
obtained by Esteban et al. I|2UU4|I . and ~0.2 dex lower 
than the NLTE abundances derived by Cunha & Lambert 
lITTMIl . 

The gas+dust corrected oxygen abundances estimated by 
Esteban et al. H1H9 8 2004) using the Cunha & Lambert 
stellar abundances in the Orion association, seem to 
be too high compared with our derived abundances 
(although still marginally consistent within the uncer- 
tainties). This result suggests a lower dust depletion 
factor of oxygen than previous estimations for the Orion 
nebula. A revision of the silicon, magnesium and iron 
stellar abundances in the Trapezium cluster stars is 
presently under way in our group for confirming this 
result. 
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